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Low lattice thermal conductivity of 
stanene 
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A fundamental understanding of phonon transport in stanene is crucial to predict the thermal 
performance in potential stanene-based devices. By combining first-principle calculation and phonon 
Boltzmann transport equation, we obtain the lattice thermal conductivity of stanene. A much lower 


: thermal conductivity (11.6 W/mK) is observed in stanene, which indicates higher thermoelectric 


efficiency over other 2D materials. The contributions of acoustic and optical phonons to the lattice 


: thermal conductivity are evaluated. Detailed analysis of phase space for three-phonon processes shows 
: that phonon scattering channels LA + LA/TA/ZA +> TA/ZA are restricted, leading to the dominant 


contributions of high-group-velocity LA phonons to the thermal conductivity. The size dependence 
of thermal conductivity is investigated as well for the purpose of the design of thermoelectric 
nanostructures. 


Research in 2D materials has experienced an explosion of interest due to their potential for integration into 
next-generation electronic and energy conversion devices! *. Recently, a new 2D material, stanene, has been 
realized on the Bi,Te; (111) substrate by molecular beam epitaxy?, and attracted tremendous interest due to its 
extraordinary properties. Stanene, a monolayer of tin film with buckled honeycomb structure like silicene, is 
predicted to be an example ofa topological insulator, which can support a large-gap 2D quantum spin Hall (QSH) 
state ?. It has been reported that a nontrivial bulk gap of ~0.1 eV opens when the spin-orbital coupling (SOC) 
is present in stanene, and the edge states are gapless with bands dispersing inside the bulk gap and helical with 
spin-momentum locking, which can be used as dissipationless conducting "wires" for electronic circuits^?. The 
free standing stanene with ~0.1 eV SOC gap is a promising applicant for room temperature and even high temper- 
ature application. In addition, a recent study has found that the thermoelectric (TE) figure of merit zT is strongly 
size dependent in stanene, and can be improved to be significantly larger than 1 by optimizing the geometric size, 


: which provides great opportunities in the next-generation TE devices”. 


All these electrical and thermoelectrical applications of stanene are closely related to its thermal transport 
properties. For instance, high TE performance requires the system to be a bad conductor for phonons but a good 
conductor for electrons, which often conflict with each other! P. In 2D QSH systems such as stanene, the edge 
states form perfectly conducting channels for electrons, leading to great improvement in the TE figure of merit 
zT!?/54, Moreover, zT of stanene can be further improved by optimizing the geometry size to maximize the 
contribution of the gapless edge states, which significantly helps to optimize the electron transport. However, 
the phonon transport properties of stanene have not been investigated in detail. Furthermore, although the size 
dependence of zT and the Seebeck coefficient in stanene has been studied for TE applications, the size dependence 
oflattice thermal conductivity, which is also crucial when designing TE nanostructures, remains uninvestigated. 

Generally both electrons and phonons contribute to the total thermal conductivity. In this paper only the 


: lattice thermal conductivity of stanene is under consideration. We calculate the lattice thermal conductivity « of 


stanene using first-principle calculation and an iterative solution ofthe Boltzmann transport equation (BTE) for 
phononsP. A much lower thermal conductivity (11.6 W/mK) is observed in stanene compared to other 2D mate- 
rials, which makes stanene a promising candidate to realize high thermoelectric performance. We also extract the 
contribution of each phonon mode to investigate the underlying mechanism behind the novel thermal transport 
properties of stanene, and find that LA phonons contribute most to the thermal conductivity (75796) over a 


: wide temperature range. It can be attributed to their high phonon group velocities and restricted phase space 


for three-phonon processes. Finally, we obtain the representative mean free path (MFP) of stanene that is useful 
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Figure 1. (a) Top view and (b) side view of stanene. 
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Figure2. Phonon dispersion and DOS of stanene. Three acoustic and three optical phonon modes are 
indicated with different colors. 


when designing TE nanostructures. The MFP in the small-grain limit is also obtained when the nanostructuring 
induced phonon scattering dominates. The role of boundary scattering in stanene nanowires is examined as well. 


Results and Discussions 

Crystal structure and phonon dispersion. Figure 1 shows the optimized structure of stanene. A similar 
low-buckled configuration is also found in silicene, germanene and blue phosphorene®*!%"’, which is different 
from the planar geometry of graphene!. The optimized lattice constant and buckling height are 4.67 A and 0.85 A, 
respectively, which are in good agreement with the previous work$7?. There are two tin atoms in one unit cell, 
corresponding to three acoustic and three optical phonon branches. 

Figure 2 shows the calculated phonon dispersion. Similar to a-Sn and 8-Sn”, the longitudinal acoustic (LA), 
transverse acoustic (TA), and z-direction acoustic (ZA) branches of stanene are linear near the I point, which 
is distinctly different from the quadratic ZA-branch dispersion in graphene”. This difference between graphene 
and stanene comes from structural difference. In planar graphene, the quadratic ZA-branch dispersion near the 
T point can be explained as the consequence of the complete decoupling of the Z modes and XY modes due to 
the Dg, point-group symmetry”). However, the buckled structure of stanene leads to the point group symmetry 
of Dza which breaks the reflection symmetry (z — —z), and therefore the Z modes of stanene couple with the XY 
modes”, which results in the linear dispersion of ZA phonons near the I point. 

Using the phonon spectrum, we calculate the group velocities of acoustic phonons along T-M and T -K respec- 
tively, as shown in Fig. 3(a,b). The anisotropy in group velocities along the two directions is weaker than that of 
black phosphorene and SnSe?*^5, As shown in Fig. 3(c), the phonon group velocities indicate weak anisotropy 
within the entire first Brillouin zone (BZ). Much smaller sound velocities are found in stanene (1300-3600 m/s 
compared to 5,400-8,800 m/s in silicene?$, 4200-6,800 m/s in MoS; 7", 4,000-8,000 m/s in blue phosphorene", 
and 3,700-6,000 m/s in graphene??). For LA, TA and ZA phonons, the sound velocities at long-wavelength limit 
are about 3548.7 m/s, 1811.4 m/s and 1303.7-1536.8 m/s, respectively. 

The Grüneisen parameter, which provides information on the anharmonic interactions, is given by 


m“ 9wjq) 
?;(q) (a) Ja (1) 


where dy is the equilibrium lattice constant, j is phonon branch index, and q is wave vector. Figure 4 shows the 
Griineisen parameters within the first BZ. Strong anisotropy is found in the Griineisen parameters of TA phon- 


ons at high frequencies, while the isotropy in the Grüneisen parameters of LA and TA phonons at low frequen- 
cies indicates that the anharmonicity of these phonons is similar from one phonon wave vector to another. The 
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Figure 3. Phonon group velocities along the (a) T-M and (b) T -K direction, and (c) within the first BZ. 
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Figure 4. Calculated Griineisen parameters within the first BZ. 


low-frequency phonons contribute most to the lattice thermal conductivity of stanene due to their higher group 
velocities and weaker scattering. Therefore, the orientation-independent phonon group velocities and Grüneisen 
parameters at low frequencies are crucial to understand the phonon transport in stanene. 


Lattice thermal conductivity. The intrinsic lattice thermal conductivity & of naturally occurring stanene 
is calculated using both the iterative and single-mode relaxation time approximation (SMA) solution, as shown 
in Fig. 5. In contrast to the anisotropy of thermal conductivity observed in black phosphorene and SnSe?*2»75, 
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Figure 5. Calculated lattice thermal conductivity « of stanene as a function of temperature ranging from 
100 K to 800 K using the iterative (red) and SMA (blue) solution of the BTE. The inverse relationship of « and 
temperature (& ~ 1/T) is displayed as dot lines. 


which possess the hinge-like structure, the isotropic & is found in stanene, which is similar to other 2D hexagonal 
materials with C; symmetry such as silicene, blue phosphorene, graphene, and monolayer MoS, 7304 

The thermal conductivity from the iterative solution at 300 K is 11.6 W/mK and 11.9 W/mK for naturally 
occurring and isotopically pure stanene respectively. Naturally occurring Sn consists of 0.97% "Sn, 0.65% ™*Sn, 
0.36% !^Sn, 14.7% !!6Sn, 7.7% !7Sn, 24.3% !!3Sn, 8.6% !?Sn, 32.4% !?°Sn, 4.6% !?Sn, and 5.6% !74Sn. 

The measured room-temperature & of 500-nm-thick tin film is 46 + 4.2 W/mK while for the thickness of 
100nm, & is 36 + 2.88 W/mK”, which are lower than the room-temperature « of bulk tin, 64 W/mK. Note that 
when the thickness of the Sn thin film decreases, phonons become more sensitive to boundary scattering’, which 
further reduces the value of «. Thus a smaller calculated & for stanene is a reasonable expectation. The Debye 
temperature Op can be calculated from the highest frequency of normal mode vibration (Debye frequency) 
Vin = 1.51THz, 


Op = hv,,/kp (2) 


where h is the Planck constant, and kz is the Boltzmann constant. The calculated Debye temperature is 72.5 K, 
which is nearly three times lower than that of 3-Sn of 200 + 3 K**. Since the Debye temperature is a measure of the 
temperature above which all modes begin to be excited and below which modes begin to be frozen out?*?, this 
low O, in stanene also indicates that the « of stanene is much lower than that of bulk tin. 

According to the kinetic formula for the thermal conductivity of an assembly of total specific heat C, average 
phonon velocity v, and MFP J, i.e. k= 1/3Cvl, when the temperature is well above the Debye temperature Op, the 
C is constant, while the dominant three-phonon anharmonic processes at high temperatures lead to lœ 1/9599, 
We fit the calculated & above © p, which reveals an inverse relation to temperature, i.e. K~ 1/T. As shown in 
Fig. 5, the fitted « (plotted as dot lines) coincides perfectly with the calculated x. Such inverse relationship of 
& and T is also found in other materials such as graphene, black phosphorene, PbSe, and Mg, (Si, Sn) when 
T» Oj5^*** Furthermore, we compare the thermal conductivity calculated from the iterative and SMA solu- 
tion. The difference between the two approaches is very small. It should be noticed that the SMA assumes that 
individual phonons are excited independently, which has no memory of the initial phonon distribution, therefore 
the SMA is inadequate to describe the momentum-conserving character of the Normal processes and it works 
well only if the Umklapp processes dominate’. Our results indicate the dominant role of the Umklapp processes 
in phonon-phonon interactions over a wide temperature range. 

Compared to other 2D hexagonal materials, the & of stanene is smaller than that of silicene (26 W/mK)”, 
monolayer MoS, (34.5 + 4 W/mK)**4 and blue phosphorene (78 W/mK)”, while at least two orders of mag- 
nitude smaller than that of graphene (3000 W/mK)???!, The lower « of stanene is due to lower Debye temper- 
ature (72.5 K compared to 500 K for monolayer MoS, and blue phosphorene!^*, and 2,300 K for graphene**). 
The Debye temperature reflects the magnitude of sound velocity. Lower Debye temperature results in decreased 
phonon velocities as mentioned above, and higher phonon scattering rates since more phonon modes are active 
at a given temperature???^, 

High TE performance requires the system to be a bad conductor for phonons and a good conductor for 
electrons. As a 2D topological insulator with lower lattice thermal conductivity, stanene can realize much higher 
thermoelectric efficiency than other 2D materials, which makes it a promising candidate for next-generation 
thermoelectric devices. 

Furthermore, we examine the contributions of different phonon branches to « of stanene, as shown in Fig. 6. 
The contribution of LA phonons is more than 57% at temperatures ranging from 100K to 800K, which is larger 
than 28% in MoS, *”, 26% in blue phosphorene" and 8% in graphene“. The contribution of ZA phonons to the 
k of stanene (13%) is much smaller than that of graphene (8096)??, blue phosphorene (44%)!”, and monolayer 
MoS, (39%)*?, but larger than silicene (7.5%)°°. It has been reported that the large contribution of the ZA mode 
to the « of graphene is due to a symmetry selection rule in one-atom-thick materials, which strongly restricts 
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Figure 6. Percentage contribution of each phonon mode to « as a function of temperature. 


anharmonic phonon-phonon scattering of the ZA mode“, while the buckled structure of stanene, silicene, and 
blue phosphorene breaks out the out-of-plane symmetry, in which the selection rule does not apply. As a result, 
the predicted contribution of ZA phonons to « in stanene is much smaller than that in graphene at 300K. To 
provide more physical insight, we investigate the phonon scattering mechanism in detail. 


Relaxation time and phase space for three-phonon processes. Considering the significant differ- 
ence in the contribution of acoustic phonon modes to the « in stanene and other 2D hexagonal materials, it is 
worthwhile to investigate the relaxation time of each phonon mode as a function of frequency, as shown in Fig. 7. 
In Fig. 7(a), the three-phonon relaxation time of LA mode is the longest, and the relaxation time of the optical 
phonons is much shorter than the acoustic phonons. We also calculate the isotopic relaxation time, as shown in 
Fig. 7(b). Isotopic scattering processes are wavelength-dependent. The long-wavelength phonons can transport 
all the heat with very little isotopic scattering”, thus we observe the relative long isotopic relaxation time of acous- 
tic phonons at low frequencies. The inverse of total relaxation time is a sum of contributions from anharmonic 
three-phonon scattering and isotopic scattering???5, 


1/7;(q) = 1/7,(q)" + 1/7,(a)”. (3) 


Thus the three-phonon processes with much shorter relaxation time will dominate the total relaxation time. 
The total relaxation time of stanene is shown in Fig. 7(c). It is in the same order of magnitude as that of graphene, 
MoS, and silicene?54959. As shown in Fig. 7(c), the total relaxation time of LA mode has large values near the 
acoustic-optical frequency gap, while it is still comparable to that of other acoustic phonon modes at low fre- 
quencies. To investigate the phonon scattering mechanism more deeply, we calculate the allowed phase space for 
three-phonon processes P3. 

In the temperature range where three-phonon processes are dominant, the total phase space for three-phonon 
processes P; is defined by?! 


(4) 


where Q is a normalization factor, and 


Pe) =D f dqD\* (q) 
J 


and 


D™ (q) = 2 dq'ê(w;(q) + wy(a’) — vj,(q + q" — G)) (6) 
JJ 


where py (q) corresponds to absorption processes, i.e. w;(q) + w(q’) = wj,(q + q' — G), whereas pr! (q) 
corresponds to emission processes, i.e. w(q) = w(q') + wj,(q — q' — G). According to Eq. (4-6), P, contains 
a large amount of scattering events that satisfy the conservation conditions and can be used to assess quantita- 
tively the number of scattering channels available for each phonon mode. Less restricted phase space for 
three-phonon processes implies a larger number of available scattering channels. Consequently, there is an inverse 
relationship between P; and the intrinsic lattice thermal conductivity of a materiaD9?!, 
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Figure 7. The (a) three-phonon, (b) isotopic and (c) total relaxation time of each phonon branch as a function 


of frequency. 
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Figure 8. Frequency-dependence of three-phonon scattering phase space at 300 K for (a) absorption processes, 
(b) emission processes and (c) all three-phonon processes. 


The P; for absorption and emission processes are shown in Fig. 8(a,b), respectively. The phase space for 
absorption processes p(? of LA phonons is much smaller than p(? of ZA and TA phonons, while the phase space 
for emission processes p(? of LA phonons is comparable to pi? of TA phonons at low frequencies but much 
larger than that of ZA phonons. According to Eq. (4), the calculated total three-phonon phase space P, for LA 
phonons is smaller than P, of ZA/TA phonons in a large part of the low-frequency range, shown in Fig. 8(c). The 
reason for the relatively small P; of LA phonons is due to the fact that the energy of LA phonons are relatively high 
and subsequently the scattering channels, LA + LA/TA/ZA — TA/ZA, are restricted. As a result, the relatively 
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Figure 9. (a) Cumulative lattice thermal conductivity of stanene as a function of the phonon MFP at 150K, 
300 K, 450 K, 600 K and 750K. The curves fitted by Eq. (7) are plotted with dot lines. (b) Thermal conductivities 
of stanene nanowires along [100] direction as a function of width. 


small P; alongside with relatively large group velocities of LA phonons as mentioned above, determines that LA 
phonons contribute dominantly to & of stanene. 
As for ZA phonons, while the relatively low energy of ZA phonons severely restricts the corresponding phase 


space for emission processes P{~’, it also provides a larger number of available scattering channels, ZA + TA/ 
ZA + LA/TA, which leads to larger pt? of ZA phonons, and hence a reduced contribution to & at low 
frequencies. 

It is well known that the dominant contribution of ZA phonons to the « of graphene is due to a symmetry 
selection rule in such one-atom-thick materials, which strongly restricts anharmonic phonon-phonon scatter- 
ing of the ZA mode“. The selection rule in graphene that arises from the reflection symmetry perpendicular 
to the graphene plane restricts the participation of odd number of ZA phonons in three-phonon processes, e.g. 
ZA + ZA © ZA, ZA + LA/TA e LA/TA. As a result, 60% of the phase space of ZA phonons is forbidden by the 
selection rule, which leads to the dominant contribution of the ZA mode to « in graphene. However, for the case 
of buckled stanene considered here, the reflection symmetry is broken, which means the selection rule does not 
apply. Consequently, the relatively large P, of ZA phonons leads to less contribution from ZA phonons to « of 
stanene, i.e. 1396 in stanene, compared to 8096 in graphene at 300K. 


Size dependence of &. Since the thermoelectric figure of merit zT can be improved by optimizing the 
geometry size to maximize the contribution ofthe gapless edge states, it is also essential to investigate the phonon 
transport in stanene nanostructures and achieve the optimized transport properties by reducing thermal conduc- 
tivity through nanostructuring in the practical design of TE nanodevices. We investigate the size dependence of & 
by calculating the cumulative thermal conductivity with respect to the maximum MFP allowed. The cumulative 
thermal conductivity of stanene at different temperatures are plotted in Fig. 9(a). The cumulated « keeps increas- 
ing as MFP increases, until reaching the thermodynamic limit above a length Lap which represents the longest 
mean free path of the heat carriers?^??, The Lapat 150 K, 300K, 450K, 600K and 750 K are 107.2 jum, 50.9 um, 
35.1 um, 25.5 um and 24.2 um, respectively. It is found that phonons with MFPs below 73.9 um at 150K, 16.9 um 
at 300K, 9.5 jum at 450K, 8.4 jum at 600 K and 6.6 jum at 750K, respectively, contribute around 75% of the total x, 
which indicates that the design of nanostructures can be utilized to reduce « and enhance zT. In order to properly 
reduce « for thermoelectric applications, we can fit the cumulative « to a uniparametric function? 


k(l < l Kmax 


se) = , 


lo 
au n (7) 


where &,,,, is the ultimate cumulative thermal conductivity, lmay is the maximal MFP concerned, and |, is the 
parameter to be evaluated. The fitted curves at different temperatures are plotted in Fig. 9(a), reproducing both 
the slope and the position of the calculated data well when Lmax > lo. It yields the parameter lọ of 5.72 um, 2.71 jum, 
1.77 um, 1.31 um and 0.95 um for 150K, 300 K, 450 K, 600 K and 750 K, respectively, which can be interpreted as 
representative of the MFP of relevant heat-carrying phonons in stanene. 

Furthermore, when the size gets smaller than l, the nanostructuring-induced phonon scattering becomes 
dominant over the three-phonon scattering, and the small-grain-limit reduced & becomes proportional to a con- 
stant value lsc". We calculate the ratio of the thermal conductivity to the thermal conductivity per unit of MFP 
in the small-grain limit; for stanene the Isg is found to be 195.0 nm, 92.1 nm, 60.9 nm, 45.6 nm and 36.5 nm at 
150 K, 300 K, 450 K, 600 K and 750K, respectively. This quantity is crucial to the thermal design for modulating 
the thermal conductivity in the small-grain limit, for example nanowires. 

It is worthy to mention that for a nanowire system, phonons with long MFPs will be strongly scattered by the 
boundary, thus the contribution from which to & will be limited. As shown in Fig. 9(b), « decreases with decreas- 
ing width of stanene nanowires, and drops to half the maximum « in the thermodynamic limit at widths about 
1.3 um. 


In summary, we predict the lattice thermal conductivity of stanene using first principle calculations and an 
iterative solution of the BTE. A much lower thermal conductivity is observed in stanene than other 2D materials. 
Our results indicate that stanene, as a 2D topological insulator with low thermal conductivity, can realize much 
higher thermoelectric efficiency and is a promising candidate for next-generation thermoelectric devices. The 
predicted contribution of LA phonons to thermal conductivity in stanene is larger than 57% due to high phonon 
group velocities and restricted phase space for three-phonon processes. The representative MFP of stanene is 
obtained for the purpose of future design in TE-application nanostructures. The MFP in the small-grain limit is 
also obtained for designing TE devices when the size decreases and the nanostructuring-induced phonon scat- 
tering dominates. 


Methods 

All the calculations are performed using the Vienna ab-initio simulation package (VASP) based on density func- 
tional theory (DFT)**. We choose the generalized gradient approximation (GGA) in the Perdew-Burke-Ernzerhof 
(PBE) parametrization for the exchange-correlation functional. We use the projector-augmented-wave potential 
with 4d electrons of tin described as valence, and a plane-wave basis set is employed with kinetic energy cutoff 
of 500eV. A 15 x 15 x 1k-mesh is used during structural relaxation for the unit cell until the energy differences 
are converged within 10-5eV, with a Hellman-Feynman force convergence threshold of 10-*eV/A. To eliminate 
interactions between the adjacent supercells a minimum of 15 A vacuum spacing is kept. 

In the calculation of phonon dispersion, the harmonic interatomic force constants (IFCs) are obtained using 
density functional perturbation theory (DFPT), which calculates the dynamical matrix through the linear 
response of electron density?*. The 5 x 5 x 1 supercell with 5 x 5 x 1k-mesh is used to ensure the convergence. 
The phonon dispersion is obtained using the Phonopy code with the harmonic IFCs as input®. 

To calculate thermal conductivity of stanene, anharmonic third order IFCs are needed, besides the harmonic 
IFCs mentioned above. The same 5 x 5 x 1 supercell and 5 x 5 x 1 k-mesh are used to obtain the anharmonic 
IFCs, and an interaction range of 8.5 A is considered herein, which includes fifth-nearest-neighbor atoms. We 
calculate the thermal conductivity of naturally occuring and isotopically pure stanene by solving the phonon 
Boltzmann transport equation (BTE) using the ShengBTE code!**®*’, which is completely parameter-free and 
based only on the information of the chemical structure. 
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